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ABSTRACT 


Evidence is presented that indicates that Monterey Sub- 
Marine Canyon was once the terminus of a major land drainage 
system. This pre-existing drainage system is not in evidence 
today because it has been altered by displacement along the 
San Andreas Fault. A numerical model based on conservation 
of mass and plate tectonic reconstructions is utilized to 
reconstruct the topography of the region as it appeared prior 
to onset of motion along the San Andreas Fault. Model results 
indicate that the Colorado River may have drained into 


Monterey Bay during early Miocene time. 
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Monterey Submarine Canyon is the largest of the California 
Submarine canyons. Traced seaward from its head near Elkhorn 
Slough in Monterey Bay, the canyon axis follows a meandering 
path for approximately 100 km before emerging as a fan valley 
Peemomacpthn Of about 3 km (Figure 1). From this point, the 
fan valley can be traced for another 400 km across the Monterey 
deep-sea fan. The volumes of the canyon and the fan have been 


3 eng. oul, O00 km> respectively (Menard, 


estimated to be 450 km 
1960). The walls of the canyon are approximately 1,500 m high 
at one point, producing relief comparable to that of the Grand 
Beamon Of the Colorado River (Figure 2). 

Such is the magnitude of this feature that two of its 
tributaries have been named. Carmel Canyon branches south 
and heads very close to shore off Monastery Beach in Carmel 
Bay. Soquel Canyon branches northward near the canyon head 
,and extends toward the town of Soquel, east of Santa Cruz. 
The major canyon axis follows a fault contact near the 
head and is flanked by Tertiary sediments on both sides. At 
its outer limits; beyond about 20 km from its head at Elkhorn 
Slough, there is no indication of faulting along the canyon 
axis and the canyon course is purely erosional (Greene, 1977). 
Soquel Canyon is cut entirely in poorly indurated seimentary 
strata, while Carmel Canyon is eroded in granodiorite on the 


east and highly indurated sedimentary rocks on the west. 
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The Monterey and Carmel Submarine Canyons 
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(Diagram 


Cpemene scent nalsGcallrornia Coast. 
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Figure 2. 


onterey 
Canyon 


Grand Canyon 
of the Coiorado 


A comparison of the profiles of the Grand 
Canyon of the Colorado River and the Monterey 
Submarine Canyon, showing them to have similar 
relief. Elevations relative to sea level are 
given in feet. Vertical exaggeration is 5X. 


(After F. P. Shepard, Submarine Geology, 3rd 
Eee an ees ana Row, 1973.) 





Like many of the submarine canyons along the California 
coast, Monterey Canyon does not lie seaward of a large land 
drainage system. The Salinas River is the largest of 
Several small rivers which drain into Monterey Bay, all of 
which appear diminutive in comparison with Monterey Canyon. 

An explanation of this enigma was proposed by Martin and Emery 
(1967) when they suggested that Monterey Canyon had received 
drainage from the Great Valley of California, through the San 
Francisco Bay region during late Miocene, Pliocene and early 
Pleistocene time. Although some questions were answered by 
this hypothesis, others remain unanswered. Martin and Emery 
noted that the amount of material in Monterey Fan could not 
have resulted from the Great Valley connection which they 
proposed. In addition, they noted the existence of buried 
erosional features which led them to conclude that Monterey 
Canyon 1S a re-excavation of a pre-existing submarine canyon. 
As Martin and Emery noted, these buried erosional features, 
Elkhorn Erosion Surface and Pajaro Gorge, pre-dated their 

Great Valley Connection. The age and size of this pre-existing 
_ canyon were established by Greene (1977) using seismic profiling 
techniques. Greene's profiles confirm the fact that this pre- 
existing canyon existed prior to early Miocene time and 1s 
still not entirely re-excavated in some areas (Figure 3). 
Greene's evidence indicates that the lower reaches of Monterey 
Canyon have been displaced northward and are represented today 


by Pioneer and Ascension Canyons. According to Greene's theory, 
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displacements along the Palo Colorado-San Gregorio Fault and 
the Ascension Fault over the past 20 million years may have 
moved these canyons into their present positions. 

Although it is presently one of the largest submarine 
Canyons in the world, there is evidence to indicate that 
Monterey Canyon is an incomplete re-excavation of a larger 
Canyon. The mechanism involved in excavation of this earlier 
Canyon and creation of Monterey Fan is unknown. 

One possible explanation is that Monterey Canyon was the 
terminus of a large land drainage system which existed in 
pre-early Miocene time. Large scale deformation of eee 
Since early Miocene has removed the canyon from its source 
and produced changes in the drainage patterns of the area 
(Clark and Rietman, 1973; Greene, 1977; Martin and Emery, 
1967; Starke and Howard, 1968). 

The topographic changes which have occurred since early 
Miocene along the west coast of North America have been exten- 
Sive. About 30 million years ago, during Oligocene time, the 
East Pacific Rise came in contact with the North American 
_ Plate. Relative motion between the Pacific and North Ameri- 
can plates came to be expressed along right lateral strike- 
slip faults on the continental margin. During early Miocene 
time, the zone of strike-slip faulting between the Pacific 
and North American Plates shifted inland to the San Andreas 
Fault. Since its inception, approximately 300 km of right 


slip has occurred along the San Andreas Fault (Crowell, 1962). 
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Such large-scale motion would undoubtedly have had great im- 
pact on any pre-existing drainage system. Reconstruction of 
topography as it existed prior to the onset of motion on the 
San Andreas Fault is difficult. Motion along the fault has 
not been constant in either speed or azimuth over this inter- 
val. In addition, compressional forces along the fault have 
uplifted areas resulting in the creation of new topographic 
features (Greene, 1977). 

Due to the complexity of the problem, a numerical model 
1s utilized to reconstruct the topography of the area as it 
appeared in early Miocene time. The model simply reverses 
motion along the San Andreas Fault utilizing data from plate 
tectonic reconstructions and imposing mass conservation. It 
does not reverse the effects of erosion or deposition, motion 
along other faults, or any other forces. However, since motion 
along the San Andreas Fault has contributed greatly to the 
alteration of topography since the early Miocene, results from 
such a model may be expected to provide useful insights into 


early Miocene topography. 


jae 
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The model was designed to simulate a reversal of the 
motion which has occurred along the San Andreas Fault over 
the past 21 million years. In simulating this reversal, 


mass must be conserved. The conservation of mass states 


that: 
dm. om, s, Vm + my - v = O (5) 
at dt : 
where 
m = mass, 
t = time, and 
vy = velocity. 


In other words, the total change in mass is equal to the local 
time change plus mass flux divergence. Since the total change 
is zero, mass 1S conserved and the local change must balance 


EranSport : 


Since total mass equals mass density (09) times total volume 


(V), then: 


aoe ee af : ¥ (pV) Ps ovV Vv (3) 
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If we assume that density is constant in space and time, 


it can be eliminated: 
aV = = > . > 7 > . _— 
Je = Vv VV VV V (4) 


Finally, when dealing with a unit area, volume is ex- 


pressed in height above/below sea level (h) (Equation (5). 


es eS he (5) 


QJ 


From Equation (5), with an initial height field and a 
velocity field, topographic changes in time may be predicted. 
The height field is a 70 x 64 array of topographic and 

bathymetric heights as illustrated in Figure 4. The grid 
aiterval is 10 minutes of longitude in the x direction and 
Mimmrnutes of latitude in the y direction (15.2 and 18.5 km, 
respectively). This interval in the x direction corresponds 
to true grid spacing at 35 degrees north, and introduces a 
departure from a spherical earth of approximately 73 at the 
. northern and southern extremes. 

While the height field simply provides initialization for 
the model, the velocity field defines the fault and drives 
the model. As illustrated in the initial velocity field 
(Figure 5), the portion of the grid which lies to the east 
of the San Andreas Fault is held fixed while that portion to 


the west of the fault is moved southeastward. Along the fault 
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itself, there 1s a gradation of the velocity field as depicted 
maeet@oure 8. At the position of the fault, the magnitude of 
the velocity field is reduced by one-half. At the next grid 
point to the east, the velocity is one-quarter as large as 
Memtsoein the moving block. At the first grid point to the 
west of the fault, velocity is three-quarters as large as 

the velocity of the moving block. This configuration results 
in a zone of velocity gradation across the fault and prevents 
the numerical instabilities which may arise from modelling 
the fault as a more severe velocity discontinuity. Previous 
attempts without a gradation of velocity across the fault 
resulted in unreasonable topography. The components of the 
velocity field, speed and direction, will be discussed 
separately. 

Directions were derived from an algorithm which generates 
smooth fields to match directions at three points where direc- 
tions are known from plate tectonic reconstructions (Silver 
fats, Upub. data). Directions of relative motion between 
the North American and Pacific Plates relative to three points 
- On the North American Plate during three time intervals are 
presented in Table I. Corresponding values from the model 
algorithm are displayed in parentheses for comparison. Since 
the model simulates a reversal of motion, all directions are 
meaueed by 180° to produce the velocity fields depicted in 
Figures 5-7. These figures illustrate the velocity fields 


used at different times in the model. The vectors to the left 
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TAB io 
Pao Wmouacemcre ah, -9 78) 


Azimuth of Average Movement Direction of Pacific 
Plate Relative to Points on North American Plate 


Values from model algorithms in parenthesis 


Age (my) Boon i I Sow Se ye o 2 W Ace Namen 2 - W 

Be-4.5 Ble 32 ieee) Poo oh.) 93h ie7? (311.77) 
ei — 1.0 Bee (328.0°) See 25.07) 819257 }§(319.5") 
mve=2).2 Ber 0 (339.2%) See est. 2°) 3235.97 (324.07) 


Side of the figures represent the moving block. It can be 
Seen that there is a very slight change in the direction from 
one end of the field to the other in all three figures as 
would be expected from Table I. On the right side of each 
figure is an area with no vectors. This is the stationary 
portion of the grid. The shorter vectors at the intersection 
of the stationary block and the moving block represents the 
fault zone. As previously discussed, there is a gradient in 
the magnitudes of the velocities across the fault. This is 
, reflected in the comparatively shorter vectors in the fault 
zone. The actual position of the San Andreas Fault 1S sSuper- 
imposed on these figures for comparison. It can be seen from 
Figures 5-7 and from the values in Table I that the motion 
changes with time to a more westerly direction becoming less 
compressional and more nearly aligned with the fault axis. 
Speeds were derived from a combination of sources. Plate 


tectonic reconstruction (Atwater, 1973) produces an accelerating 
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velocity field which would have resulted in approximately 

610 km of offset if all motion had been expressed as fault 
displacement. How much of this displacement is actually 
expressed in offsets along the San Andreas and other faults 

is unclear. Evidence for a 260 km offset along the San Andreas 
Fault (Crowell, 1962) establishes a conservative estimate 

and is the reference displacement utilized in this model. Al- 
though displacements from plate tectonic reconstruction do not 
match displacements derived from geologic evidence, the accelera- 
tion of motion between two plates, as derived from plate tec- 
tonic reconstruction, is not disputed by geologic evidence 
(eewater, 1973). This acceleration of motion 1s incorporated 
into the model by reducing Vanloen ties from plate tectonic 
reconstruction such that the total resultant offset equals 260 
km. This is accomplished by multiplying the velocities by 

a scaling factor, s = 260 km/610 km = .426. Resultant scaled 
velocities as well as original velocities derived from plate 


tectonic reconstruction (Atwater, 1973) are listed in Table II. 


TABLE II 
Time Period Relative Motion Relative Motion 
(Million years Pacific/N. Ameri- San Andreas 
before present) can Plates Fault 
21--10 ier See CIly wae OMG Sem 35 
10--4.5 4.0 cm/yr iL ietnvalea 
4.5-- en OE eV ata g 
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The height and velocity fields are substituted into the 
Semeriuity Equation (Equation (5)) using finite differencing 
schemes. In discussing differencing schemes, the following 
coordinate system and symbology will be used: 

a =r Lerenco.mouwnt, Neight at a location on 


aa ©) a horizontal plane identified by coor- 
dinates i (east) and j (north) at time t 


u = velocity component in the i (east) 
direction 
v = velocity component in the j (north) 
daeeci ion 
At = interval of time (200,000 years in this 
model) 
Ax = interval of space on the i axis. For 
example, the distance between H,. . 
alg ol ash oar Seba Joe) 
ere de) 
Ay = interval of space on the j axis. The 
distance between Hai, d,t) and Hi, je+l,t) 


The primary centered differencing scheme used in this 


model expresses Equation (5) in the following way: 


et), €+1) ie eG aates L) = Gee ll ine) H(i, 9,6) ) 
: Zit 7 Die 


eevee im ol Cea jl t) | 
ZA; 


Peer meme Clee ise), t) | 
ag id 


= eWay jtl,t)—-v(i,j-1,t) ] 
2ZAy 


(6) 





This scheme is centered in space and time. It can not be 
used initially since it requires height fields at two differ- 
ent times (t and t-l1) in order to predict the height at ttl. 
In addition, it can't be used on the grid boundaries since 
it would require the input of values beyond the range of the 
data field. However, it offers the advantage of centered 
differencing and favorable numerical stability. 

Forward differencing must be used to predict the height 
fields on the first iteration for subsequent input into the 
primary centered scheme. Using forward time and centered 


space differencing, 


Peete, ttl) -H(1,7j,t) om oe OR GUA aS), PEG leet ee) 
At _ Py Re 


a vee fre) elec, pee eee (ale ie Sa 
2Ay 


Peco sina st) uate) 
ZK 


Pee eee piney ea Vil tot, ©) | 
2Ay 


(7) 


As with all numerical methods, numerical stability is an 
important consideration and limits the size of the time step 
eeeyeto be utilized. The Courant-Friedrichs-Levy (CFL) 
mematerom (Courant et al., 1928) for computational stability 


(Equation (8)) applies for this model. 





Wien ec 
| nee | < a (8) 
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It controls the size of the time step (At) and grid 
spacing (Ax,Ay) such that motion does not cover more than one 
grid space between computations. With this model, At must be 
less than or equal to 500,000 years. The time step used is 
200,000 years. 

Boundary conditions are illustrated in Figure 4. As pre- 
viously discussed, centered space differencing can not be 
used on the perimeters of the grid. Two different types of 
boundary conditions are used in this model. Boundaries in the 
direction of motion of the moving block (the eastern and southern 
boundaries) must be permeable. For this reason, the Upstream 
Differencing scheme (Haltiner and Williams, 1980, p. 130) is 
used on these boundaries. This scheme produces a permeable 
boundary where it intersects the moving block and establishes 


meee 2a DOuNndary condition for the stationary block: 


Seummemee (i, j,t) = .  uli-l,j,t) [H(i,j,t)-H(i-1l,j,t) ] 
ive i Ax 


= Pe een, ce) aa) a,c) | 
Ax 


(8) 


Fixed boundary conditions are imposed on the northern 
and western boundaries. Maintaining the points on these 
boundaries at their initial heights introduces errors where 
these points lie on the moving portion of the grid. As 
topography on these boundaries moves toward the interior 


of the grid, a false field is created in their previous 
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positions. Resultant erroneous values are eliminated from 
output by moving the northern and western boundaries inward 
(for display only) such that the northwest corner of the grid 
Maintains its position relative to the moving block. This 

is illustrated in Figure 9. 

The model produces output at 3 million year intervals. 
Output is displayed in a three dimensional form as well as a 
contoured plan view at each interval (Figures 10-25). The 
three dimensional computer graphics simulate viewing the ter- 
rain from a vantage point at high elevation over the Pacific 
Ocean looking northeast. Submarine bathymetry is included and 


all vertical heights are exaggerated by a factor of 15. 
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THE NORTHERN AND WESTERN 


BOUNDARIES MOVE INWARD 
SUCH THAT THE NORTHWEST 
CORNER OF THE GRID MAIN- 
TAINS ITS POSITION REL- 
ATIVE TO A POINT (R) ON 
THE MOVING BLOCK, 





Figure 9. Movement of boundaries for display 
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Figure 10. Topography 21 million years before present 
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Paogure 12. Topography 18 million years before present 
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Figure 14. Topography 15 million years before present 
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Figure 16. Topography 12 million years before present 
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Figure 18. Topography 9 million years before present 
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Figure 22. Topography 3 million years before present 
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In viewing the figures depicting model reconstructions 
of paleotopography, the design of the model must be kept in 
mind. Coarse resolution, as well as the simplifying assump- 
tions of the model limit the scope of consideration to large 
scale effects. 

Figures 10-25 are model outputs generated at three 
million year intervals. At each interval there is a three- 
ditiensional mepogeaniic display and a contour map, both 
constructed from model results for that particular time. 

In all three-dimensional plots, the coastline (based on 
present sea level) is represented by a heavy solid line. 
The points corresponding to the cities of Monterey, Bakers- 
field, and Los Angeles have been annotated for reference. 

All contour maps are oriented with north at the top. 
Peeenymetry appears as dashed contours while all contours >0 
appear as solid lines. The San Andreas Fault 1s represented 
as a heavy solid line. Inferred drainage patterns are 
“indicated by a heavy dashed line. The points corresponding 
to the cities of Monterey, Bakersfield, and Los Angeles are 
labelled M, B, and L respectively. 

Although the model runs backward in time, this discus- 
Eon will start with model output at 21 million years before 


present and proceed foward in time for simplicity. As 
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Suteput for each time interval is discussed, reference will 
= made to both figures (3-D and contour) representing 
mfewut for that interval. 

As expected, topography at 21 million years before 
present (Figures 10 and 11) differs greatly from present- 
day topography. Monterey Bay lies far to the south-east 
of its present position, approximately 100 km west of 
Bakersfield. Bakersfield, and points to the east of the 
San Andreas Fault, appear much as they are today. The Great 
Valley of California extends northwest of Bakersfield and is 
bounded on the east by the Sierra Nevada Mountains. The 
Mojave Desert lies to the south of the Sierra Nevadas. 
Along the fault itself, topography is depressed, particularly 
in the area to the southwest of the Mojave Desert. Areas to 
the west of the fault appear radically different than they 
do today. All topographic features which correlate to today's 
topography appear far to the southeast of their present posi- 
tions. Although the general shape of the coastline and the 
mountain ranges are recognizable, they appear different than 
‘they do today. Monterey and San Francisco Bays are different 
in shape and location but are easily recognizable. The Coast 
Ranges to the north appear at lower elevations than they do 
today but are not greatly different in the area to the south 
of Monterev. The Transverse Ranges are somewhat lower and 
they are located to the southeast of their present positions 


as are the Peninsular Ranges. The position of the Peninsular 
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Ranges correlates geographically with the present-day 
Meettion of the Gulf of California. 

Despite the fact that this reconstruction is the product 
of a relatively simple geological model which does not take 
into account the effects of erosion, deposition, or displace- 
ments along other faults, the large scale features which 
appear in these figures are supported by geologic evidence. 
Seemeereetra Nevada Mountains (Curtis et al., 1958) and the 


Colorado Plateau (Eardley, 1962) are features which came 


Gnto existence long before Miocene time. Also, the Gulf of 
Seecooniad did not exist at this time (Larson et al., 1968; 
Larson, 1971; Van Andel and Shor, 1964). The depressed 


areas along the fault southwest of the Mojave Desert corre- 
late geographically with the Salinas, Caliente, San Joacquin, 
Ridge, and Soledad baSins. The geologic history of these 
ffetms (Norris and Webb, 1976; Blake et al., 1978) also 

shows good chronological correlation with the model 
meeonstruction. 

It is apparent from Figures 10 and 11 that drainage 
‘patterns were much different in early Miocene time. Drainage 
from the Colorado Plateau (Colorado River) could not flow 
Beuthiward into the Gulf of California as it does today. 
Rather this drainage flowed westward and entered the ocean 
somewhere to the north of the Transverse Ranges. Exactly 


where the terminus of this drainage system was located can 
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not be determined from model results alone. However, evidence 
previously presented indicates that Monterey Bay was the 
terminus of a major land drainage system at this time (Greene, 
meee The inferred drainage is illustrated in Figure 11. 

How long this drainage pattern had been in existence prior 

to early Miocene time is a question beyond the scope of this 
study. However, insights into changes which have occurred 
Since early Miocene may be obtained by examining model output 
at intervals over the last 21 million years. 

As expected from the relatively low velocity of motion 
along the fault from 21 million years before present to 10 
million years before present (Table II), there is little 
change in topography over this period (Figures 10-17). It 
appears that the drainage pattern previously established 
would have been preserved over this interval. Although there 
is evidence to indicate that Monterey Submarine Canyon was 
filled and re-excavated during this time (Greene, 1977), this 
could have been the result of sea level fluctuations and does 
not necessarily imply large scale changes in land drainage 
‘patterns. 

From 10-4.5 million years before present, motion along 
the fault increased in speed and became slightly less com- 
pressional in azimuth (Tables I and II, Figure 6). Topography 
at 9 million years before present (Figures 18 and 19) shows 


that compression along the fault in the area southwest of the 


48 





Mojave Desert is closing the depressions in that area. At 
the same time, the Peninsular Ranges have been moving north- 
westward, leaving a depressed area in the southeastern corner 
Sethe Grid. At this point it is not possible to determine 
drainage patterns, as illustrated in Figure 19. At 6 million 
years before present (Figures 20 and 21), it is clear that 
compressional motion along the fault has uplifted the area 

to the southwest of the Mojave Desert to such an extent that 
this area no longer serves as a conduit for drainage. Simul- 
taneously, the depressed area in the southeast corner has 
expanded due to continued movement of the Peninsular Ranges 
out of this area. It 1S not possible to determine that 
drainage actually flowed southward at this point since 
topography beyond the southeastern corner of the grid is not 
presented. However, geologic evidence (Larson et al., 1968; 
Larson, 1971) indicates that the Gulf of California developed 
at about this time. For this reason, it is assumed that 
drainage from the Colorado Plateau flowed southward much as 
it does today. 

Velocity of relative motion is slightly higher over the 
last 4.5 million years to the present (Table II) and there 
Mena slight change in the azimuth of direction (Table I and 
Figure 5). However, the trends concerning drainage continue 
during this interval as depicted in Figures 22-25. North- 


westward displacement of points to the west of the San 
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Andreas Fault, compression in the area of the Transverse 
Ranges, and expansion of the Gulf of California continue to 


the present and are continuing today. 





iv. CONCLUSIONS 


In arriving at conclusions based on the results of a 
numerical model, care must be taken to consider the inherent 
weaknesses of the model and to examine available evidence 
from other sources for verification. In presenting model 
results, only large scale features and general trends were 
discussed. It is clear, from the description of the model, 
that these are the limits of its credibility. However, 
within these limits, model results correlate well with 
available evidence and indicate that the Colorado River was 
the erosional force involved in the excavation of Monterey 


eanyon. 
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